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Abstract 

We examine the Inert Doublet Model in light of the discovery of a Higgs-like boson with a 
mass of roughly 126 GeV at the LHC. We evaluate one-loop corrections to the scalar masses and 
perform a numerical solution of the one-loop renormalization group equations. Demanding vacuum 
stability, perturbativity, and 5-matrix unitarity, we compute the scale up to which the model can 
be extrapolated. From this we derive constraints on the model parameters in the presence of a 
126 GeV Higgs boson. We perform an improved calculation of the dark matter relic density with 
the Higgs mass fixed to the measured value, taking into account the effects of three- and four- 
body final states resulting from off-shell production of gauge bosons in dark matter annihilation. 
Issues related to direct detection of dark matter are discussed, in particular the role of hadronic 
uncertainties. The predictions for the interesting decay mode h° — >> 77 are presented for scenarios 
which fulfill all model constraints, and we discuss how a potential enhancement of this rate from 
the charged inert scalar is related to the properties of dark matter in this model. We also apply 
LHC limits on Higgs boson decays to invisible final states, which provide additional constraints on 
the mass of the dark matter candidate. Finally, we propose three benchmark points that capture 
different aspects of the relevant phenomenology. 
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I. INTRODUCTION 



The Inert Doublet Model (IDM) is the simplest among the models with two Higgs dou- 
blets. In addition to the Standard Model (SM) particle content, it contains an extra doublet 
of complex scalar fields which couples to the SM scalar and gauge boson sector but not 
to the fermions. Moreover, it involves a discrete Z 2 symmetry under which the new scalar 
doublet is odd and all the other particles are even, which makes that the new "inert" doublet 
particles can only appear in even number in interaction vertices. 

The IDM was first introduced more than three decades ago in studies of electroweak 
symmetry breaking (EWSB) [1]. Long after, it was proposed as a model that can provide a 
viable dark matter candidate according to the thermal relic picture @, El , since the neutral 
scalars contained in the new doublet can be seen as weakly interacting massive particles 
(WIMPs) and play the role of the dark matter (DM) in our universe. Due to its rich phe- 
nomenology for cosmology and particle physics, the IDM has received considerable attention 
jH]. Its DM candidate captures all the basic mechanisms through which the observed relic 
density can be generated in WIMP models Q: The "correct" relic abundance [6] can be 
achieved by adjusting couplings, by approaching or taking distance from resonances, or by 
co-annihilating with another particlejj 

Additionally, the IDM was advocated to allow for a heavier SM-like Higgs boson com- 
patible with electroweak precision tests, m h o > 200 GeV, without resorting to unnatural 
fine-tuning [3] (for recent work also considering this possibility, see, e.g^ Ref. Q). The new 

|, and they provide 
911 . While by now a 
111, the model still 
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states predicted by the IDM have been subjected to collider bounds 
an interesting phenomenology for the Large Hadron Collider (LHC) J 
heavy Higgs boson with SM couplings is experimentally ruled out |10|, 
remains attractive due to its dark matter features. Furthermore, the IDM provides an inter- 
esting example of interplay between DM and Higgs physics, since the SM-like Higgs boson 
is one of the basic means of communication between the a dark"0 and "visible" sectors of the 
model. Hence, both the mass and the couplings of the SM-like Higgs boson are of crucial 
importance to assess whether the IDM can indeed explain the dark matter in the universe. 

With the recent announcement of the observation of a Higgs-like resonance with a mass 
of around 125-126 GeV by the ATLAS and CMS collaborations 0, HJJ, as well as the 



supporting hints from the D0 and CDF experiments [HI, it appears likely that the particle 
responsible for EWSB (or at least one of them) has been discovered. In this spirit, we find it 
timely and interesting to examine implications of this observation for the IDM. Apart from 
the obvious consequence that the number of free model parameters is reduced by one, other 
interesting features appear, as we shall describe in the present paper. Moreover, despite the 
attention that the IDM has received in the community, most studies rely on lowest order 
predictions only. An exception is the study beyond leading order that was performed in 
Ref. [13]. This work is, however, limited to the interesting case of radiative electroweak 
symmetry breaking a la Coleman- Weinberg, a scenario leading to rather extreme parameter 
values. Another exception is the very recent paper [HI, where higher-order corrections to 
DM direct detection in the IDM are calculated. 

In the present work, we perform an analysis of the IDM parameter space assuming that 
the LHC is indeed observing a (SM-like) Higgs boson with a mass M h o ~ 126 GeV. In 
Sec. (Til we introduce the model and the relevant notation. In Sec. [Till we present one-loop 

1 Coannihilation is absent, e.g., in the simpler singlet scalar model. 

2 This term is used in a slightly abusive way here, since the model does not contain what is usually dubbed 
a "dark sector" in the litterature. 
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corrections to the scalar masses in the IDM, and expressions for the one-loop renormalization 
group equations (RGEs) for the model's quartic couplings, which are used to study vacuum 
stability, perturbativity, and unitarity constraints. Experimental constraints from collider, 
low-energy, and cosmological data are presented in Sec. [IV] The corresponding numerical 
analysis, which contains the main results of this work, is presented in Sec. |VJ where we 
perform extensive scans over the model parameter space. The implications for dark matter 
and the interesting Higgs decay into two photons are discussed in detail. Based on our 
analysis, we identify benchmark scenarios capturing general features of the parameter space. 
Finally, Sec. IVII contains a summary of our results and the main conclusions. 



II. THE INERT DOUBLET MODEL AT TREE LEVEL 

The inert doublet model (IDM) contains, in addition to the Standard Model (SM) particle 
content, a second complex scalar doublet. The model Lagrangian is constructed so as to 
satisfy an exact Z 2 symmetry (parity) under which all SM particles, including one of the 
scalar doublets, are even and the second scalar doublet is odd. Since this discrete symmetry 
prevents mixing between the scalars, one of the doublets, denoted by i/, is similar to the 
SM Higgs doublet. Given that the second doublet, which we denote by $, has odd Z 2 
parity it is inert in the sense that its component fields do not couple singly to SM particles. 
The requirement of renormalizability then also forbids all tree-level couplings to the fermion 
sector. 

Imposing the Z 2 symmetry has two further important consequences. First, the inert 
doublet $ does not acquire a vacuum expectation value. Second, it forbids several of the 
terms appearing in the general two- Higgs- doublet model scalar potential 0. More precisely, 
the tree-level scalar potential of the IDM takes the form 

Vb = /^l^ + z^l^ + Ailtf^ (2.1) 

2 - 

In the general case, the are complex parameters. Although considering this possibility 
can have interesting consequences for CP- violation and electroweak baryogenesis [16| , in this 
work we limit ourselves to the case of real values. Upon electroweak symmetry breaking, 
the two doublets can be expanded in components as 

where v = V2 (0\H\0) « 246 GeV denotes the vacuum expectation value of the neutral 
component of the doublet H. The h° state corresponds to the physical SM-like Higgs-boson, 
whereas G° and G ± are the Goldstone bosons. The inert sector consists of a neutral CP-even 
scalar ii°, a pseudo-scalar A , and a pair of charged scalars H ± . 

A phenomenologically important consequence of the Z 2 symmetry is that the lightest 
Z 2 -odd particle (LOP) is stable. If further the LOP is either H° or A° : this (neutral) state 
can play the role of the DM candidate, analogously to the Lightest Supersymmetric Particle 
(LSP) in supersymmetric models with imparity conservation. 

At the tree level, the scalar masses are obtained from the potential in Eq. ( 12.11) . When 
the potential is expanded in the component fields, the masses of the physical states are given 
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by 

m\ =nj + 3\ lV 2 , (2.3) 

m 2 H0 =/4 + A L v 2 , (2.4) 

m 2 A0 = f4 + A<?v 2 , (2.5) 

m 2 H ± = l4 + \^v 2 . (2.6) 

Here we have introduced the useful abbreviations 

\l = \ (A 3 + A 4 + A 5 ) , (2.7) 

\s = \ (A 3 + A 4 - Ag) . (2.8) 

As in the Standard Model, we have the scalar potential minimization relation m 2 h0 = — 2fi 2 = 
2X 1 v 2 that can be used to eliminate the parameter fi\ after electroweak symmetry breaking. 
The IDM scalar sector can hence be specified by a total of six parameters 

{Ai, A 2 , A 3 , A 4 , A 5 , /x 2 }, (2.9) 

which can be exchanged through the above equations in favour of the physically more mean- 
ingful set 

{m h o, m H o, m A o, m H ±, A L , A 2 } , (2.10) 

that is often used in phenomenological applications. It is worth noting that the parameter 
Xl has a simple physical interpretation as the h° — H° — H° coupling at the tree level 
(similarly, A# is the h° — A — A coupling). In contrast, the parameter A 2 , which is common 
to both parameter sets, appears only in quartic self couplings among inert particles and does 
therefore not enter any physically observable process at the tree level. However, it plays a 
role once the theory is considered beyond leading order. In the following, we shall make use 
of both sets of parameters given here. 



III. THE INERT DOUBLET MODEL BEYOND THE TREE LEVEL 

Most of the work that has been performed within the IDM so far, such as collider sig- 
nals, vacuum stability considerations or relic abundance constraints, has been based on 
pure leading-order calculations. Generally speaking, it is meaningful to go beyond tree level 
when theoretical uncertainties of some calculation start to become comparable to the corre- 
sponding experimental uncertainties. In many models this is the case, for example, for the 
well-measured relic abundance of dark matter in our universe. Another reason can be to 
investigate effects which are absent at tree level but arise in leading order at the one-loop 
level as it is the case, e.g., for the decay h° 77. 

In the IDM, radiative corrections to inert scalar annihilations, which are relevant for 



DM phenomenology, are expected to be rather small 17]. Here, strong couplings do not 
intervene in the relevant Feynman diagrams (with the exception of a final state vertex 
correction involving gluon exchange among final state quarks), contrary to the Minimal 
Supersymmetric Standard Model, where 0(a s ) corrections can have a sizeable impact in 



this context [18|. Moreover, the number of diagrams with electroweak couplings at one-loop 
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order is rather limited, again unlike the case of the MSSM, where electroweak corrections 
can become significant due to the very large number of contributing diagrams [191 ]. 

With the electroweak scale being probed by both colliders and dark matter detection 
experiments, the different new physics models are constrained and the allowed magnitude of 
their couplings becomes limited. In particular, the IDM parameter space is known to present 
three mass regimes for the lightest odd particle (LOP), where the observed relic density can 
be obtained 0, H, (^[H]: the so-called low- (1 GeV < M LO p < 10 GeV), intermediate- 
(10 GeV < M LO p < 150 GeV), and high-mass (M LO p > 500 GeV) regimes (the distinction 
being somewhat arbitrary). As we shall see in the following, the low- and intermediate- 
mass regimes are already quasi-excluded from the XENON100 experiment results [23]. As 
a consequence, in order for the IDM to still provide DM- compliant scenarios in this mass 
range, one has to rely on resonances and thresholds to achieve the correct relic density 
without violating the null XENON results, which would happen for sufficiently large coupling 
values. In particular, we will see that the only points of the interesting intermediate mass 
regime that survive the constraint are hang close to the Higgs "funnel" region and the 
Vl^Vl^* production threshold (see also 0, 3] for recent works on the subject), i.e. the region 
where LOP annihilation is dominated by nearly on-shell h° exchange or gauge boson pair- 
production. 

In this regime, given the accuracy in the experimental determination of the relic density, 
a mass shift of the order of 2 GeV is sufficient to render a particular point in parameter 
space excluded or notJl If radiative corrections are sufficiently large to shift the LOP mass 
for which the correct relic density is observed in or out of the viable region, then the com- 
putation of higher-order corrections can become important. This is well-known in the case 
of supersymmetric models such as the MSSM or the NMSSM, see, e.g., Ref. [24j|. The exact 
values of the involved masses are also known to be essential in the case of co-annihilation 
of the LOP with the next-to-lightest Z 2 -odd particle (NLOP). In this case, a shift in one or 
both masses can change the importance of this particular channel and thus the prediction 
of the relic abundance. 

In order to perform arbitrary calculations at the one-loop level within the IDM, we have 
implemented the particle content and the interactions (in 't Hooft-Feynman gauge) in model 
files for the Mathematica package FeynArts/FormCalc fiBj]. The implementation has been 
performed in parallel and cross-checked with the help of the two packages FeynRules f26j 
and SARAH Q. 



A. One-loop corrections to the scalar masses 



In the Feynman-diagrammatic (FD) approach, the one-loop corrected scalar masses are 
obtained by evaluating the corrections to the two-point functions. Denoting the mass for a 
generic scalar field (j) appearing in the renormalized Lagrangian by m^, the corresponding 
(renormalized) inverse propagator, T^, can be written as 

r (p 2 ) = i [ P 2 - K - ReS *(p 2 ))] , (3-1) 

where S</>(f> 2 ) is the (again, renormalized) self-energy of the field (j). From this equation, two 
further concepts can be defined. The first is that of the running mass, which we denote 

3 Note that this mass difference is also similar to the current difference between the values for the Higgs 
mass extracted by ATLAS and CMS in different channels. 
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V f 
FIG. 1. Generic classes of Feynman diagrams contributing to the one-loop self-energy of scalars in 
the Inert Doublet Model. The upper diagrams involving scalars and vector bosons (V) contribute 
to the masses of both the SM-like (H) and the inert scalar doublet (<£>), while the lower diagrams 
involving a vector boson (V) or a fermion (/) are only present for the SM-like Higgs doublet. 



by m^(Q). As indicated by the argument, this mass carries an explicit dependence on the 
renormalization scale Q. The running mass is obtained from the inverse two-point function 
evaluated at zero external momentuno 



m$(Q) = ilW = 0) = m\ - Re E (p J = 0). (3.2) 

The second type of one-loop mass is obtained by solving Eq. ( 13.11) for the value of p 2 = M| 
that gives T(M|) = 0. In practice, this equation is solved iteratively with usually very fast 
convergence. Since it corresponds to a pole of the propagator, is refered to as the pole 
mass. Which one-loop mass definition should be used depends on what is most convenient 
in a given situation. Working in fixed order of the perturbative expansion, the impact of 
this choice on physical observables is a higher-order effect. Below we shall make use of both 
the running mass and the pole mass when studying different aspects of the IDM. 

The one-loop self-energies appearing in Eq. (13.11) generally receive contributions from 
scalars, gauge bosons, and fermions. In the IDM, the diagrams involving fermions are 
restricted to the SM-like Higgs boson, which is also the only state that couples singly to the 
vector bosons. Generic classes of diagrams contributing to the scalar self-energies are shown 
in Fig. [TJ Working in the MS renormalization scheme, we evaluate these in the 't Hooft- 
Feynman gauge, which leads to the expressions for the running masses given in Appendix 
lAl Before showing numerical results, we should also mention that these results have been 
checked against those presented in Ref. [13], and overall agreement was found (modulo the 
fact that Ref. [l3| does not include the 0(a) corrections coming from vector boson loops). 



Moreover, the SM-type contributions for m^o have been checked against [28| . 

We find that the relative difference between the MS and the tree-level masses can, de- 
pending on the exact parameter values and the chosen renormalization scale, reach up to 
30% for the SM-like Higgs boson h°, and up to around 10% for the inert scalars. In Fig. 
El this is illustrated for an example scenario U As expected, the difference between the pole 
and the running MS masses is smaller. For the example in Fig. El it amounts to up to 2 GeV 
(about 1.5%) for h° and up to almost 1 GeV (about 1.5%) for H°, which is representative 



4 An alternative method which yields the same expression for the running masses is through the (one-loop) 
effective potential, see [13] for the case of the IDM. 

5 This scenario corresponds to the benchmark point I defined in Tab. HT1 below. 
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FIG. 2. Running one-loop masses in the MS scheme of the SM Higgs boson (h°) and the lightest 
inert scalar (LOP) as a function of the renormalization scale Q (blue, solid lines). This example 
corresponds to the input parameters m h o = 125.3 GeV, m H o = 66.0 GeV, m^o = 131.0 GeV, 
m H ± = 137.0 GeV, A^ = 0.12, and A2 = 8.6 • 10 -4 . In the upper panels, the running mass is 
obtained by interpreting the input mass as the tree-level mass (green, dash-dotted), while in the 
lower panel the input value is interpreted as the pole mass (red, dashed). In this example, the 
LOP is the CP-even scalar H°. 



for typical parameter points in the IDM. The situation is similar for the masses of the other 
inert scalars, A and H ± ) which are not shown here. 



B. Renormalization group equations for the quart ic couplings 

Renormalization group equations (RGEs), when evaluated at one-loop accuracy, resum 
leading scale logarithms to all orders. They are therefore useful for studying the evolution 
of parameters under renormalization over a large energy range. To treat all parameters at 
a similar footing, we find it more convenient here to work in the basis of quartic couplings, 
Eq. ( 12.91) . rather than the mass basis. The renormalization group equations for the quartic 
couplings (i = 1 ... 5) are driven by the corresponding beta functions, f3\ i: according to 
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The latter receive contributions from scalars (s), gauge bosons (g), and fermions (y). The 
different pieces of j3\ i are obtained by diagrammatic calculation^ The results for the different 
contributions are given in Appendix [Bl Our expressions are in agreement with earlier results 



for the general 2HDM [30, [3jJ (recalling that in the IDM A 6 = A 7 = 0). 



IV. CONSTRAINTS 



The various constraints which can be imposed on the (tree-level) parameters of the inert 
doublet model (IDM) have been extensively discussed in the literature. A summary is, e.g., 



given in [32|]. For completeness, we dedicate this Section to discuss which constraints are 



relevant here and how they are treated in our analysis. 
A. Theoretical constraints 

A first class of constraints on the inert doublet model (IDM) comes from the minimal 
requirement for a physically sensible and (perturbatively) calculable theory. First, it is 
required that the couplings should not be larger than some value which makes a perturbative 
treatment meaningless. We thus demand all quart ic scalar, gauge, and Yukawa couplings to 
fulfill a common constraint 

M, |^|, \ Vi \ < K. (4.1) 

A minimal condition is that this inequality be respected at the input scale where the pa- 
rameters on the left-hand side are specified. The value of K which should be imposed for 
a valid perturbative expansion is then somewhat process-dependent. Considering the more 
stringent requirement that Eq. (14.11) remains valid under renormalization group evolution, 
we can use Eq. (13.31) for guidance. Since the coefficients appearing in the beta functions 
for the scalar quartic couplings are all 0(1), it is clear that an instability (Landau pole) 
will appear for |A^(Q)| > An in the absence of accidental cancellations. In practice, smaller 
values of |A^| can also lead to instabilities, but we choose to be conservative and impose 
Eq. (13.31) with K = An as the perturbativity limit at all scales. 

A second requirement is that the scalar potential, Eq. (12.11) . should be bounded from 
below and that upon electroweak symmetry breaking (EWSB) it develops a minimum which 
renders the electroweak vacuum stableQ For this requirement to hold we demand [H, 34 

Ai(Q),A 2 (Q) >0, 

A 3 (Q) > -2VAi(Q)A 2 (Q), (4.2) 

a 3 (Q) + HQ) - MQ)\ > -2 V / A 1 (g)A 2 (g) 

at the scale Q. 

A final constraint comes from the requirement that the scattering matrix (5-matrix) of 
every quantum field theory must be unitary. In the case of weakly coupled theories, it 
is sensible to require that the tree-level scattering matrix elements satisfy unitarity limits, 

6 The renormalization group equations can also be computed by demanding scale invariance of the one-loop 
effective potential. For an example of this approach for the singlet scalar model, see Ref. [291 ]. 

7 The criterion dicussed here implies absolute stability of the symmetry-breaking vacuum. More relaxed 
constraints would be obtained by also allowing for metastable configurations with a sufficiently long 
lifetime, see Ref. [33I for a recent discussion for the SM. 
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which corresponds to imposing upper bounds on them. For the general 2HDM, the bounds 
were first derived in |35|. Here, we use the form for the eigenvalues of the scalar and vector 
scattering matrices of [36], and we require them to be smaller than 167T. This corresponds 
to saturation of the unitarity limit with the tree-level contribution. 



B. Oblique parameters 

In models where the dominant effects of new physics appear as corrections to self energies 
of the (SM) gauge bosons, the effects can be parametrized in terms of the three "oblique" 
(Peskin-Takeuchi) parameters S, T, and U [37], which vanish in a pure SM calculation. In 
the 2HDM, the contributions to the U parameter are negligible, which makes it convenient 
to work in the approximation U = 0. This assumption has been verified explicitly in our 
numerical analysis. For the case of U = 0, recent experimental limits on the remaining two 
parameters are [38] 

S = 0.06 ± 0.09, T = 0.10 ± 0.08. (4.3) 

These values are based on a reference (SM) Higgs mass of m™* = 120 GeV and a reference 
top mass of m\ ei = 173 GeV. We impose the limits resulting from Eq. (14.31) at the 2 a 
confidence level as constraints on the IDM contribution. 



C. Collider searches 



The first constraint from direct searches at colliders comes from the invisible decay width 
of the Z boson. If the decay mode Z —> H°A° is open, the subsequent decay A H°ff 
(or H° —> A°ff for the inverse mass hierarchy) would lead to Z decay events with fermion- 
antifermion pairs (//) and missing energy in the final state. A detailed analysis has shown 
that this decay is incompatible with LEP data, which implies that the decay width of 
Z — > H°A° must be small. It is convenient to implement this constraint as {3, 0] 

M H o+M A o > M z . (4.4) 

Assuming a fixed mass hierarchy M H o < M^o , a more detailed analysis of the IDM parameter 
space with respect to LEP data leads to the limit M^o > 100 GeV [8]. Considering both 
possible mass hierarchies between H° and A , we require 

max{M if o,M A o} > 100 GeV. (4.5) 

Finally, limits on the mass of the charged scalar can be obtained by considering their po- 
tential pair production and subsequent decay into neutral Higgs bosons at LEP. Converting 
existing limits on the search for charginos and neutralinos, which present the same final state 
topology at colliders, leads to the bound M H ± > 70 — 90 GeV [39|]. For practical reasons, 
we adopt the intermediate limit 

M H ± > M w . (4.6) 
In order to have a neutral DM candidate, we include as a final requirement 

M H ± > M LOP = min(M H o, M A o) (4.7) 

for defining viable points in parameter space throughout our analysis. 
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D. Dark matter relic density 



The recent results from the WMAP satellite, combined with other cosmological measure- 
ments, constrain the dark matter relic density to ^cdm^ 2 = 0.1126±0.0036 [6]. In most of 
the numerical analysis below, we require the IDM relic density (considering either H° and 
A as the dark matter candidate) to respect this limit within 3 a, i.e. we demand Q^oph 2 to 
lie within the interval 

0.1018 < Q LOF h 2 < 0.1234. (4.8) 
In some cases we shall however relax this requirement and only impose the upper bound, 

^Loph 2 < 0.1234, (4.9) 

which corresponds to the situation where the IDM contribution is only partly responsible 
for the observed DM density. 

A subtle point in the calculation of the relic density concerns the intermediate mass 
regime, and in particular the mass interval 50-80 GeV, where the IDM is known to produce 



values for ^ L op^ 2 in the correct range. It has been pointed out [2l|, [22| that in this mass 
regime, which lies close to the WW final state threshold, contributions to the total self-/co- 
annihilation cross section coming from three-body final states can be substantial or even 
dominant. In particular the Vl^Vl^* contribution (and to a lesser extent also ZZ*) should be 
taken into account to obtain reliable predictions. For this purpose we use a modified version 
of micrOMEGAs [40] which includes three-body final states' contributions. This micrOMEGAs 
version also accounts for four-body final states with two virtual gauge bosons. We include 
these contributions in our analysis, and we find that they are not always negligible; the 
effect on the computed relic density can be up to a few percent H For more details on the 
implementation of the relic density calculation we refer to the corresponding manual |4l| . 



E. Dark matter direct detection 

Further constraints stem from (null) direct searches for dark matter. Until recently, the 
most stringent bounds on the spin-independent WIMP-nucleon scattering cross-section came 
from the XENON10 jll] and XENON100 jll] experiments for the low- and intermediate- 
mass regimes respectively. Recently, the XENON100 experiment presented updated con- 
straints on the spin-independent cross-section, which strengthen the previous bounds by a 
factor of roughly 2-6 |23| . We impose these more stringent bounds in the numerical analysis 
that follows. 

Our calculation of the relevant spin-independent WIMP-nucleon cross sections is again 
performed using the micrOMEGAs package. It is well-known that predictions of direct de- 
tection rates are not without uncertainties. One among the most important sources of 
uncertainty (at least from the particle physics side) is related to the strange quark form fac- 
tor of the nucleon, f Ts: which loosely speaking describes the "strange quark content" of the 
nucleon. This quantity is usually parametrized in terms of the so-called pion-nucleon sigma 
term g^n (3). The value of g^n relies on experimental data that is at present poorly known. 
The values which are widely used come mostly from lattice QCD and chiral perturbation 



theory calculations |44|, [45 



8 We have noted small differences with respect to the results presented in Ref. [21|, which we attribute to 
the inclusion of the doubly virtual final states. 
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Hadronic form factors, and notably fas, are of particular importance when the dominant 
mechanism for WIMP-nucleon scattering is Higgs exchange (as is the case in the IDM), 
since the Higgs boson couples preferentially to heavy quarks and the strange-quark density 
is sufficiently large to provide a non-negligible contribution to the scattering cross section. 
Since the results of varying f Ts can be quite dramatic (13, Si, 

we perform our calculations 



adopting two distinct values for this parameter: f Ts = 0.2594, corresponding to the com- 
monly used value g^n = 55 MeV, and fr s = 0.014, a value obtained recently from first 
principles through lattice QCD methods [47|. It should be noted that, as clearly stated 
in Ref. [43], this study does not treat systematic uncertainties. It nonetheless provides us 
with an indicative value for f Ts: englobing numerous recent results that point to much lower 
values than previously estimated. 



V. NUMERICAL ANALYSIS 
A. Setup and strategy 

To analyze which regions of the IDM parameter space are compatible with theoretical 
and experimental constraints, we perform a numerical analysis scanning randomly over the 
model parameters. Since the first step of our analysis will involve RGE running of the 
parameters starting from the input scale (which we take to be the pole mass M z of the 
Z-boson), we define the input in the MS scheme in the basis of the Lagrangian parameters 
{A 2 , A 3 , A4, A 5 , with the exception of the Higgs boson pole mass M h o. For the IDM- 
specific parameters, we consider the following input ranges: 

< A 2 (M Z ) < 2, 
-2 < A 3A5 (M Z ) < 2, (5.1) 
-v 2 < i4(M z ) < (ITeV) 2 . 

We select the distribution of ji\ in order to obtain a uniform distribution of the CP-even 
scalar mass fn H o. Moreover, in order to improve the sampling efficiency, we generate more 
points in the regions of parameter space where points are expected to remain valid up to 
high scales. The chosen ranges of Eq. (15.11) are found to yield both a good parameter space 
coverage and a sufficiently efficient scan with an acceptable number of viable points, both 
from the point of view of theoretical constraints at the input scale and with respect to the 
experimental constraints that we shall impose in the following. 

For the first part of the analysis (for which the results are shown in Fig. [3j), we do not 
make any assumption on the value of the Higgs boson mass and let it vary freely within 
the bounds < M h o < 500 GeV. We remind the reader that here, and in the following, 
capital letters always refer to one-loop pole masses. A sample of 10 6 parameter space points 
is generated within the aforementioned ranges, with the only requirement that they fulfill 
vacuum stability, perturbativity, and unitarity at the input scale Mz- 

For the second scan, which is used throughout most of the following analysis, we make 
the assumption that the LHC experiments have indeed observed a (SM-like) Higgs boson 
with a mass around M h o ~ 126 GeV. We therefore fix the mass of the SM-like Higgs boson 
in the IDM to lie in the range 

M h o = 125.7 ±0.6 GeV, (5.2) 
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which is obtained by a Gaussian combination of the two experimental measurements jlOl . Ill 
This also fixes the value of Ai at the input scale. We then sample 10 7 points in the IDM 
parameter space, {A 2 , A 3 , A4, A 5 , again demanding that the stability, perturbativity, and 
unitarity constraints are satisfied at the scale M^. 

A further important SM input is the top-quark mass, which in both scans is sampled 
from the range [481] 

M t = 173.2 ± 0.9 GeV. (5.3) 
This mass is also interpreted as the pole mass, related to the running mass at one loop by 



m t {m t ) M t 



3^.(M«) 



(5.4) 



To solve the renormalization group equations numerically, we use a modified version of 
the 2HDMC code (version 1.2) [49]. At each scale, the conditions for perturbativity, vacuum 
stability, and unitarity are evaluated. As soon as one of the three conditions fails, we record 
the corresponding scale, as well as the failing condition. The maximal scale we consider is 
Q = 2 • 10 16 GeV (the "GUT" scale), at which the evolution is terminated and surviving 
points are also recorded. 2HDMC is also used to evaluate the oblique parameters (at the input 
scale), and we apply the bounds from direct searches at colliders as described in Sect. UVCl 
For each parameter space point, we evaluate the scalar pole masses from the running pa- 
rameters as described in Sec. IIII A[ This set of parameters is passed to micrOMEGAs, which is 
used to compute the relic density and the WIMP-nucleon spin-independent scattering cross 
section. As a final step, we compute the Higgs decay rates (again using 2HDMC) to check 
whether an enhancement of the important h° —> 77 decay mode (as hinted by the LHC) 
could be reproduced in any of the viable parameter points, and to apply constraints from 
invisible Higgs decays into LOP pairs. 



B. Extrapolation scale 

Our first goal is to examine the maximal scale A up to which the IDM can be extrapolated 
while remaining a consistent and calculable quantum field theory, in the sense of preserving 
stability of the electroweak vacuum, perturbativity of all couplings, and unitarity of the 
scattering matrix as discussed above. 

We first employ the scan with variable Higgs mass. For comparison, we also perform the 
same excercise for the SM, where the only free parameter is taken to be the mass of the Higgs 
boson h°. Both analyses are performed on the same footing, using one-loop renormalization 
group equations and varying the top quark mass within the limits of Eq. (15.31) The results 
are presented in Fig. [3J where we show the maximal extrapolation scale of the SM (left) and 
the IDM (right) as a function of M h o. The value corresponding to the LHC observation, as 
given by Eq. (15.21) . is indicated by the dashed line. 

In the left-hand side plot of Fig. [3j we can see the two well-known SM bounds coming 
from perturbativity (red), and vacuum stability (black), constraining the Higgs mass from 

9 Far more sophisticated analyses exist for the SM, where state-of-the-art calculations rely on three- loop 
RGEs for the couplings [33| . They also consider the possibility of sufficiently long-lived metastable vacua. 
This analysis shows that the SM Higgs mass range yielding an absolutely stable vacuum is enlarged by 
a few GeV compared to the lowest order prediction, while the recent LHC discovery points towards a 
(long-lived) metastable vacuum for the SM. Such a detailed analysis, while crucial, goes far beyond the 
scope of the present work. 



12 




FIG. 3. Left: Results from a one-loop RGE analysis of the maximum extrapolation scale A which 
satisfies perturbativity (red), stability (black), and unitarity (not relevant) as a function of the 
Higgs mass M h o in the Standard Model. The value corresponding to the LHC signal is indicated 
by a dashed line. Right: The corresponding results in the Inert Doublet Model for the parameter 
variations according to Eq. (15.11) . 



above and below, respectively. The width of these two lines corresponds to the uncertainty 
from variation of the top quark mass. 

The introduction of the additional inert doublet significantly modifies the extrapolation 
constraints as a function of the Higgs mass, as can be seen on the right-hand side of Fig. [31 
While the modification of the upper bound from perturbativity is less pronounced, the 
lower one from stability is greatly relaxed, with the electroweak vacuum being stable up to 
high scales for a larger choice of the higgs mass than in the SM case. This effect is due 
to the contributions of the inert scalars to the potential, that can counteract the tendency 
of the top quark Yukawa coupling to destabilize the vacuum by driving Ai negative, see 



) q i 

Refs. [50|, 51j. Interestingly, this leads to a vacuum stability constraint that only manifests 
itself fairly early in the RGE evolution. In other words, parameter space points which have a 
stable vacuum at the input scale are likely to remain stable up to high scales. It is therefore 
generally the case that one of the other two requirements (perturbativity or unitarity) is 
the more constraining for M h o (as can be seen From Fig. [3]). We find that the range for 
which the IDM Higgs mass can be extrapolated up to the GUT scale is larger than that 
in the SM, and in particular that lower values for M h o are allowed. The measured value 
M h o ^126 GeV, which at our level of approximation does not allow for the SM to be valid 
up to the GUT scale, see Fig. [3] (left), is therefore allowed in the IDM. We expect this 
conclusion to hold also in a more sophisticated analysis since, as already mentioned, the SM 
bound on M h o from stability is relaxed by the inclusion of higher-order corrections [33|. 



In the following we constrain the Higgs boson mass to the range of Eq. (15.21) suggested by 
the LHC discovery. We also impose the oblique parameter and collider constraints presented 
in Sees. IIVBI and IIVCI As before, we evolve the RGEs until one of the three theoretical 
requirements fails and record the corresponding scale. The results are presented in Fig. 
where the scale A, at which the first of the constraints fails, is presented against four different 
parameters/predictions of the model, namely M L qp, \M h o — M a q\ ) \ L) and A 2 . The colour 
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FIG. 4. Maximum scale A up to which the IDM satisfies perturbativity, vacuum stability, and 
unitarity as a function of different model parameters/predictions: the LOP mass (top left), the 
mass difference \M H o — M A o\ (top right), the H° — H° — h° coupling A^ (bottom left) and the A2 
parameter (bottom right). The criterion which failed at the lowest scale is indicated by the colour 
code: perturbativity (red points), vacuum stability (black), and unitarity (green). 



coding is the same as in Fig. O red points fail due to the perturbativity constraint, black 
points fail due to vacuum stability, and the green points are those for which unitarity is 
violated at the scale A. We observe that, as also noted previously, the vacuum stability 
constraint is relevant mostly at rather low scales A < 10 7 GeV. Unitarity violation can 
occur as late as the GUT scale. 

From Fig. [H we also learn that extrapolating the model towards the GUT scale does not 
favour a particular range for the mass of the dark matter candidate (the LOP). However, 
the coupling parameters are strongly constrained when reaching scales close to the GUT 
scale. This is due to the requirement of perturbativity, which is much easier to fulfill if the 
numerical values of the couplings are small. As a further consequence, the H° — H° — h° 
coupling X L and the four-scalar coupling A 2 have to be chosen in a rather narrow interval for 
the model to be valid up to the GUT scale. The mass difference between H° and A is also 
limited if the GUT scale validity requirement is applied. This is due to the fact that m H o and 
m^o stem from a common mass scale, /x 2 , and their difference is proportional to X^v 2 (see 
Eqs. ( 12.41) ). Since the quartic couplings are bound due to the perturbativity requirement, 
the mass difference for valid points naturally decreases when approaching the GUT scale. 
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FIG. 5. Projections of the parameter space of the Inert Doublet Model on the planes of A^ against 
A3, A4, A5, and A#. The green regions correspond to all valid points in the scan, while the black 
regions show the points which remain valid up to the GUT scale. 



Parameter 


Scan range 


Valid range 


GUT range 


A 2 


(0.0,2.0) 


(0.0,2.0) 


(0.0,0.15) 


A 3 


(-2.0,2.0) 


(-1.1,2.0) 


(-0.25,0.75) 


A 4 


(-2.0,2.0) 


(-2.0,2.0) 


(-0.80,0.90) 


A 5 


(-2.0,2.0) 


(-2.0,2.0) 


(-0.30,0.30) 


Al 




(-0.3,3.0) 


(-0.15,0.40) 


As 




(-0.5,3.0) 


(-0.15,0.40) 



TABLE I. Input ranges for the IDM quartic couplings, allowed ranges when imposing constraints 
at the input scale, and after extrapolation to the GUT scale. 

The same behaviour is observed for the remaining coupling parameters, as shown in Fig. [5] 
where the IDM parameter space is projected on the planes (A 3 , A^), (A4, A^), (A 5 , A^), and 
(As, Al). The (light) green regions in this figure correspond to parameter space points that 
pass the aforementioned constraints at the input scale M z (constraints from DM observables 
are not imposed here; they will be discussed in the next section). Thus, they constitute the 
maximal phenomenologically viable regions of the IDM within our choice for the parameter 
ranges. At the same time, we highlight in black the points corresponding to parameter 
values for which the IDM can be extrapolated up to the GUT scale. As can be seen from 
this figure, the constraints imposed at the input scale are not too prohibitive, and they are 
satisfied for essentially the full parameter ranges considered in our scan. It can also be seen 
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that the IDM can be extrapolated up to the GUT scale in sizeable regions of the parameter 
space, with the ranges of A^, As, and A 5 for which this can happen being somewhat more 
limited, whereas A3 and A4 can reach higher absolute values. In Tab. U we summarize the 
allowed ranges for the quartic couplings, both at the input scale and after the evolution to 
the GUT scale. 



C. Dark Matter 

With the previous results at hand, we now turn to the IDM dark matter phenomenology, 
assuming as previously that the Higgs and top quark mass are fixed within the ranges given 
by Eqs. (15.21) and (15.31) . respectively. In Fig. [6] we show the predicted relic density as a 
function of the DM candidate mass, M L op (left), and the triple coupling of a DM pair to 
the (SM-like) Higgs boson (A^ when the LOP is if , A# when it is A ) (right). The 3a limit 
from the 7- year WMAP data is represented by the red-dashed region. 

From the left-hand side plot in Fig. El we can see that the correct relic density can be 
achieved in the mass regimes that we described in SecJIIIl The viable parts of the low- and 
intermediate- mass regimes extend from 3 GeV up to roughly 120 GeV. From this value and 
up to approximately 500 GeV, the predicted relic density is too low and it can reach the 
WMAP levels again above 500 GeV. In the right-hand side plot, we see that the values that 

(As) can take while yielding the correct relic density lie in the range —0.4 < X L ^s J$ 0.4, 
with positive values prefered by the points viable up to the GUT scale. 

In order to better illustrate the impact of the WMAP results, in Fig. [7] we project the 
IDM parameter space onto the (M L op, A^) plane, demanding that the relic density satisfies 
only the upper 7- year WMAP limit (left), or both the upper and the lower bound (right). 
In particular the right plot shows how restrictive the full DM constraint is, and it allows to 
discuss in more detail the various mechanisms responsible for producing the correct values 
for ^lop^ 2 in the different mass regimes. 

At very low masses, below 4 GeV, LOPs annihilate dominantly into r pairs through 




FIG. 6. The dark matter relic density versus the LOP mass Mlop (left) and the coupling of a 
LOP pair to the Higgs boson A^s (right). The green points correspond to all valid points in the 
scan, while the black region shows the points which remain valid up to the GUT scale. The red 
dashed line indicates the WMAP central value of Eq. (14.8)) . 
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FIG. 7. The viable IDM parameter space projected on the (Mlop, Xl,s) plane imposing only the 
upper limit (left) and the upper and lower limits (right) of the WMAP range of Eq. (14. 8ft . The 
colour coding is the same as in Fig. [61 

s-channel Higgs exchange. The most energetic part of the LOP's Boltzmann velocity dis- 
tribution can also give some annihilation into bottom quarks. The correct relic density is 
achieved by adjusting X^s which, given the relative smallness of the r Yukawa coupling, has 
to be fixed at rather large values. These values decrease as phase space opens up, but also 
as a larger part of the WIMP Boltzmann distribution passes the bb production threshold. 
Above M L op ~ 4 GeV, the required Xl,s values drop by a factor 3-4, compensating for the 
larger bottom Yukawa coupling. Continuing up to roughly 50 GeV, the bulk of the points 
lie on a thin strip where A^s decreases mildly, corresponding to dominant DM annihilation 
into bb pairs through Higgs boson exchange. We remind the reader that due to the LEP 
direct search constraints, coannihilation is absent in this mass range. For Mlop ~ 50 GeV, 
we observe that the WMAP-allowed values for X^s decrease further, approaching X^s = 
at precisely M L op = M h o/2. Two effects are at play in this mass region. First, the LOP 
mass approaches the point where the Higgs propagator becomes resonant. For a constant 
Xl value, the LOP self-annihilation cross section then increases dramatically. As a result, in 
order to get the correct relic density, X^s must attain small values. The second important 
effect in this mass region comes from the contributions of virtual gauge boson final states, 
and especially from V^*)]/^*. Even for vanishing triple coupling of a DM pair to h° : there is 
an additional four-vertex involving two DM particles and two W bosons. This results from 
the covariant derivative terms in the Lagrangian, acting on the $ doublet. The strength 
of this vertex is thus characterized by a gauge coupling, and does not depend on any free 
parameter of the model (see also Ref. [HI]). This constitutes a distinct positive contribution 
to the DM annihilation cross section that pushes the required value for \Xl^s\ further down. 

For the IDM high-mass regime, a detailed explanation for the behaviour of the the relic 
density constraint has been presented in Ref. [20]. Above roughly 120 GeV (the precise 
value depends on M h o (22|), the cross section for self-annihilation, in particular into gauge 
boson final states, becomes prohibitively large for the lower WMAP bound to be satisfied. 
We therefore see that, whereas this region is allowed in the sense of not overshooting the 
WMAP limit, Figs. El [7| (left), it fails once we demand the IDM to fully account for the 
measured DM density of the universe, Fig. [7] (right). However, for M L op > 500 GeV, an 
interesting effect comes at play. The dominant contribution to the cross section of the type 
H°H° —> VV (A A VV) here comes from the longitudinal gauge boson components. 
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This cross section is mediated either by a direct quartic coupling, or by a i/^-channel 
exchange of H° : A , or H ± ) which scales as M^jM\ (in the case of the Z final state, and 
similarly for the W). The annihilation cross section therefore becomes very large as Mlop 
increases (with a similar dependence also on the NLOP mass). However, when the H° and 
A '/H ± are nearly mass-degenerate, there is a cancellation taking place between the t/u 
channel contributions and the four-vertex diagram. This cancellation is exact for an exact 
degeneracy. For example, with A^ = 0, M H o = 700 GeV, and M A o = 701, the relic density 
would be too high, and to satisfy the WMAP bound a non-zero value for A^ is required. In 
this manner, the WIMP depletion rate can be balanced by varying the LOP-NLOP mass 
splitting and the A^ parameter to obtain the correct mixture of transverse and longitudinal 
gauge bosons in the final state. These solutions are always found for small LOP-NLOP 
mass splittings, and require some tuning of the value of Xl,s- In practice we find that the 
maximal allowed mass splitting for the points in our scan is of the order 10 GeV. 

Further constraints from dark matter come, as we have already mentioned, from direct 
detection experiments, and most notably the latest XENON100 limits on the WIMP-nucleon 
spin-independent scattering cross section, a$i. In Fig. [8] we show the viable IDM parameter 
space in the usual (Mlop^si) space and confront the model to the latest XENON100 
exclusion bounds (23| (dashed red line). As discussed in Section HV El we adopt two distinct 
values for the strange quark nucleon form factor, and the results are shown in the left- and 
right panels of Fig. [HI respectively. In this figure all points respect both the upper and the 
lower WMAP bounds. The behaviour of the LOP-nucleon scattering cross section follows 
quite closely the corresponding behaviour of the coupling \ L ^$ in Fig. [7J since the only way 
of coupling the LOP to quarks at tree level is through i-channel Higgs exchange (the relevant 
coupling being simply Xl,s)- It can be seen clearly that the low- and intermediate- mass 
regimes are almost fully excluded. The only surviving points are those for which the correct 
relic density is achieved through a combination of small values for the A^s coupling and 
quasi-resonant annihilation to an s-channel Higgs boson and virtual gauge boson final states. 
For completeness, we note that the very low mass regime (M L op < 10 GeV) is also excluded 
by the XENON10 bounds on low-mass WIMPs [421 ]. which for the sake of simplicity has 
not been depicted here. The high mass regime, on the other hand, remains unaffected by 
current direct detection bounds. 

In Fig. [8] we can also clearly observe a consequence of fixing the Higgs mass to a constant 



FIG. 8. Projection of the viable IDM parameter space on the (Mlop^si) plane against the latest 
XENON100 limits (red dashed line) for f Ts = 0.2594 (left) and f Ts = 0.014 (right). 
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value. The resonance regime, which for a variable Higgs mass consists of an entire region, 
acquires now the well-known two-branch "funnel" structure that is present, e.g., in models 
such as the MSSM, where the Higgs mass is fixed to a very limited range. The funnel ex- 
tends quasi-symmetrically in two narrow regions around M h o/2. The LOP cannot approach 
too close to the resonance since the relic density would then be too low. With the Higgs 
mass fixed, the freedom of moving within the funnel region essentially disappears, which 
demonstrates that a greater accuracy is demanded in order to assess whether a parameter 
space point is excluded or not. A similar comment applies to the high mass regime and the 
precise tuning of the LOP-NLOP mass splitting that we discussed above. 

Moreover, we observe the anticipated, quite dramatic, impact of fys on the LOP-nucleus 
scattering cross section. As we pointed out, the driving mechanism for WIMP-nucleon 
scattering in the IDM is i-channel Higgs exchange. The Higgs couples preferentially to 
strange quarks, hence varying the strange quark form factor can have an impact of up 
to a factor 7 in the scattering cross section. While in this particular case most of the 
parameter space remains excluded, this is a transparent illustration of the importance of 
knowing f Ts with better accuracy to make precise statements. Using the more conservative 
value, we find that the mass ranges allowed by WMAP and direct exclusion bounds are 
50 GeV < M LOP < 80 GeV and M LOP > 500 GeV. 

As a final remark on direct detection, we should mention that in a recent paper the one- 
loop electroweak corrections to the LOP-nucleon scattering cross section were computed flij . 
According to the findings, these corrections give positive contributions to the scattering 
cross section of the order of 10 -11 -10 -10 pb, that can actually superseed the tree-level 
contributions in non-negligible regions of the parameter space (in particular for small 
values). These contributions do not depend on the IDM couplings, since the relevant vertices 
involve gauge couplings, and for the most should be quite insensitive to variations of /t s - The 
most striking impact is that these corrections basically set a lower value for the scattering 
cross section over the full range of LOP masses, which turns out to be within the reach of the 
upcoming XENON- IT experiment. This is a nice example of the complementarity among 
direct dark matter and collider searches, since the high mass regime, most probably out of 
the LHC reach, will be probed by the next generation of direct detection experiments. 



D. LHC Higgs phenomenology 

The phenomenology of the Higgs- like state observed at M h o ^126 GeV remains mostly 
unchanged in the IDM compared to the SM. Since the inert scalars do not have tree-level 
couplings to fermions, the SM production modes and dominant decay channels are unaf- 
fected, which leaves the rates compatible with the data from different channels measured by 
ATLAS [Tot] and CMS fill ]. On the other hand, the prediction for one of the experimentally 
most important decay modes, h° —> 77, can be affected at leading order. In the SM, the 
h° —> 77 decay is mediated by loop diagrams, where the most important contributions come 
from W bosons and top quarks, while in the IDM it receives an additional contribution 
from the charged scalarsjp Following the LHC discovery, this has also been discussed in the 
context of the 2HDM (Hlf 

We calculate the partial width T(h° 77) and the corresponding branching ratio for the 
points in the parameter space that satisfy the constraints, including the WMAP bounds on 

10 The leading-order prediction for the decay h° — )• could be similarly modified, but this mode has not 
yet been observed experimentally. 
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FIG. 9. Relative branching ratio i? 77 = BR(/i° —> 77)idm/BR(/i° —> 77)sm as a function of the 
charged scalar mass M H ± (top) and the coupling A3 (bottom), applying only the upper WMAP 
limit on the relic density (left), and both the upper and lower limits (right). The colour coding is 
the same as in Fig. El 



the relic density. In Fig. [9J we show the ratio 

BR(h° -> 77 ) IDM 
BR(h° -+ 77)sm : 



(5.5) 



which in the IDM corresponds to the modification of the inclusive h° —> 77 rate (recall that 
the h° production in, e.g., gluon fusion, remains unchanged). The colour coding used in 
Fig. [9] is the same as above, where green indicates all viable points in the scan, and the 
black points in addition remain stable and perturbative up to the GUT scale. Considering 
first the upper left plot, we show i? 77 as a function of M H ± with only the upper WMAP limit 
on the relic density applied. For large values of M H ± the rate decouples to the SM value 
(i? 77 — > 1), whereas for M H ± < 300 GeV, we see that the two photon rate can be modified 
substantially. Here we observe scenarios with an enhancement in i? 77 above the SM. The 
effect is most pronounced for low M H ± < 100 GeV, where the possible enhancement exceeds 
30%. It should be pointed out that these low values for the charged scalar mass are allowed 
in the IDM without contradicting either indirect constraints from flavour physics f53[, or the 
direct LHC searches for charged Higgs bosons [54j]; again, since the charged scalar does not 
couple to fermions. On the other hand, values for M H ± ^100 GeV would limit the mass 
for the IDM dark matter candidate, M L qp, to the low mass regime. 
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FIG. 10. Projection of the viable IDM parameter space for Mlop versus \l,s (left) and asi (right), 
with the latest XENON100 limits (red dashed line). The colour coding is the same as in Fig. El The 
full DM constraint from WMAP is applied, as well as the upper bound BR(/i° invisible) < 0.65. 



As can be seen from Fig. [9j the R 77 enhancement is not compatible with requiring a 
scenario valid up to the GUT scale, since these scenarios (black points) always lead to 
R 77 < 1. The explanation behind this is given in the lower left plot, which shows i? 77 versus 
A 3 . Here it can be seen that, once the WMAP constraints are applied, the sign of A 3 must 
be positive for scenarios to be viable up to the GUT scale. At the same time, the sign of this 
parameter governs the interference between the charged scalar contribution and the (SM) 
WW contribution, and an enhancement can only be obtained for a negative value of A3. 

The second set of plots in Fig. [9] (right) are the same as those to the left, except that both 
the upper and the lower WMAP limits on the relic density have now been applied. These 
plots show important difference compared to the results where only the upper WMAP limit 
is applied. First of all, the number of allowed points is naturally much reduced. A second, 
very important, difference is the complete absence of points with i? 77 > 1 in this case. Once 
the (full) relic density constraint is imposed, the IDM rate for h° — >• 77 should therefore be 
expected to come out similar to that in the SM. 

Another possibility, which is visible in Fig. [9J is to have more than an order of magnitude 
suppression of i? 77 . This large suppression can be caused by "invisible" Higgs decays to the 
LOP, h° H°H° (h° A A ) for M LOP < M h o/2. The resulting low values of i? 77 are not 
compatible with the observation of a SM-like Higgs boson at the LHC, which can serve as 
a further constraint on these scenarios. In addition, the ATLAS Collaboration has recently 
published upper limits (55| on the possible size of an invisible branching ratio for a SM-like 
Higgs boson. At the observed value for the Higgs mass, this corresponds to the constraint 



BR(h° -> invisible) < 0.65 (5.6) 

at 95% CL. In the IDM, when h° decays to LOP pairs is open, the invisible branching ratio 
tends to become very large. The effect of applying Eq. (15.61) as a constraint can be seen 
in Fig. [TOl which shows the valid points in the (Mlop, ^l,s) an d (Mlopj^si) planes after 
applying this bound. As illustrated by this figure, this constraint from the LHC provides 
a complementary limit to that from direct detection, in ruling out IDM dark matter in the 
low- and intermediate mass regimes. 
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M h o 


M H o 


M A o 


M H ± 


Al 


A 2 


f2 L op^ 2 


CSI 


(pb) 




A (GeV) 


I 


125.3 


66 


131 


137 


0.012 


9 • 10" 4 


0.1113 


7.1 • 


10 -io 


0.91 


> 2 • 10 16 


II 


125.8 


973 


972 


976 


0.170 


0.022 


0.1224 


4.4- 


10 -io 


1.00 


> 2 • 10 16 


III 


125.9 


91 


262 


102 


-0.359 


1.182 


7 • 10" 4 


3- 


lO" 7 


1.25 


1.2 • 10 4 



TABLE II. Scalar masses (in GeV) and quartic couplings at the input scale, Mz, for three bench- 
mark scenarios, with the corresponding values for the relic density, Qhoph 2 , direct detection cross 
section, asi, the h° — > 77 rate relative to the SM, i? 77 , and the maximum extrapolation scale, A. 



E. Benchmark scenarios 

As a final result of our analysis, we extract three benchmark points which capture different 
aspects of the IDM phenomenology. The scalar masses and quartic coupling for these points 
are summarized in Tab. [Til 

The first benchmark point has M L op = 66 GeV ~ M^o/2. It lies at the border be- 
tween the Higgs-funnel and Vl^Vl^* threshold regions, where the annihilation proceeds partly 
through a h° resonance and partly into (virtual) gauge bosons. For the given point, the com- 
bination of the two mecanisms leads to the value of ^lop^ 2 = 0.1113, which is in agreement 
with the WMAP interval of Eq. (14.81) . The second benchmark scenario, point II, lies in the 
high-mass regime, the LOP being the pseudoscalar A . Its mass is almost degenerate with 
the other inert scalars, so that the relic density of Q^oph 2 = 0.1224 is achieved through the 
interplay of pair- annihilation and co- annihilation between ^4°, H° : and H ± . Both scenarios 
I and II can be extrapolated to the GUT scale, but neither of them present an enhanced 
h° —> 77 decay rate as compared to the Standard Model. Point II is out of the reach of 
LHC, but might be probed by the XENON1T experiment. 

Finally, point III is inspired by recent LHC observations rather than by its dark matter 
phenomenology. It presents an enhanced h° —> 77 decay rate (i? 77 = 1.25), but a rather 
low relic density (fi L op^ 2 ~ 7 ' 10 -4 ). Extrapolation of this scenario is only possible up to a 
maximum scale of about 10 TeV, where the requirement of perturbativity fails. Let us note 
that in principle it would have been possible to find a scenario with i? 77 ^ 1.25 and a higher 
extrapolation scale (see e.g. Fig. [11]), but the corresponding charged scalar masses are found 
to be very close to the limit M H ± > Mw We have therefore chosen a more conservative 
point where M H ± > 100 GeV. 

The fact that it is not possible to satisfy all the constraints which have been discussed in 
our detailed analysis while obtaining an enhanced i? 77 within the IDM is again illustrated 
in Fig. [TTJ This figure shows the correlation of i? 77 to the relic density Q^o^h 2 (left) and 
the maximum extrapolation scale A (right) obtained from our parameter space scan. We 
clearly see that points with i? 77 > 1 are only found for Q L Q F h 2 < 0.1. In the high-mass 
regime, due to the requirement M H ± > M L op, the charged scalar is also rather heavy, and 
the branching ratio of the decay h° — > 77 is equal to the SM value. In the low-mass regime, 
the LOP is light enough to allow for invisible decay of the SM-like Higgs, so that its decay 
rate into photons is suppressed. In the intermediate mass region, as mentioned before, A 3 is 
restricted to positive values in order to satisfy the upper and lower WMAP limits. At the 
same time, a negative value of A 3 is necessary to achieve i? 77 > 1 (see Fig. [9]). Consequently, 
both conditions cannot be satisfied at the same time. Relaxing the lower WMAP bound 
allows for points similar to our benchmark scenario II, i.e. presenting an enhanced value of 
i? 77 > 1, but with a lower dark matter relic density ^lop^ 2 0.1. 
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FIG. 11. Correlation of the relative branching ratio i? 77 = BR(/i° — » 77)idm/BR(/i° — » 77)sm and 
the dark matter relic density ftj^oph 2 (left) and the maximum extrapolation scale A (right). The 
colour coding is the same as in Fig. El 

Similarly, the situation of i? 77 > 1 occurs only for extrapolation scales lower than A < 10 7 
GeV, as can be seen from Fig. [IT] (right). As discussed before, the extrapolating to high 
scales puts an upper limit on the absolute value of the coupling parameters. In consequence, 
the decay rate of h° — >■ 77 is smaller because of the limited quartic coupling A 3 . 

VI. SUMMARY 

We have presented an extensive study of the Inert Doublet Model in the light of the 
recent potential discovery of a Higgs-like boson with a mass around 126 GeV at the Large 
Hadron Collider (LHC). To this end, we have computed the one-loop corrections to the scalar 
masses of the model and evaluated their numerical effect. In particular, we have shown the 
difference between the pole mass and the running MS masses to be of the order of a few 
percent. This is sizeable in the context of dark matter annihilation, which often relies on 
precise mass differences with respect to a resonance, or on the presence of co-annihilation 
which is equally sensitive to the precise relations among the masses of the involved particles. 

Moreover, we examined the high-energy behaviour of the IDM when treating the Higgs 
boson mass either as a free parameter, or under the assumption that the particle observed at 
the LHC corresponds to the IDM Higgs boson. In particular, we require the perturbativity 
of all couplings, stability of the electroweak vacuum, and unitarity of the ^-matrix. By em- 
ploying the one-loop renormalization group equations, we have shown that these constraints 
become less prohibitive than in the case of the Standard Model. More precisely, the IDM 
can be extrapolated up to the GUT scale for a non-negligible range of parameter choices 
and for larger ranges of the Higgs boson mass, with the experimentally observed value lying 
firmly within the region where the IDM can remain valid up to very high energy scales. 

We have then imposed a series of constraints coming from collider (LEP, oblique param- 
eters) and dark matter (WMAP, direct detection) experiments, combining them with the 
aforementioned theoretical requirements. We find that the oblique parameter constraints has 
only mild effects upon the IDM parameter space, and the main result of the LEP bounds 
is to push the NLOP mass above 100 GeV, which in turn results in the absence of neutal 
co- annihilation in the lower LOP mass regime (roughly below 80 GeV). 
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As a further consequence, the relic density in this mass region is driven by the self- 
annihilation cross section, which can be modified by tuning the couplings of the LOP (A^ or 
Xs for H° or A , respectively) to the Higgs boson, or by approaching the Higgs resonance 
and WW production threshold. Far from the resonance/ threshold, relatively large (or 
As) values are required for the WMAP bounds to be satisfied. This also implies a large 
WIMP-nucleon scattering cross section, which brings the model in tension with the current 
XENON100 bounds. The only surviving region lies close to the Higgs resonance/ WW pro- 
duction threshold, where \l,s, an d consequently the coupling to quarks via Higgs exchange, 
has to be small in order for the depletion rate of DM not to be too high. The significant 
reduction of the available parameter space demonstrates the crucial character Higgs mea- 
surements can have for a number of dark matter models. Another measurement which we 
have used to constrain the mass region M L op < M h o/2 is the upper limit on an invisible 
branching ratio for the SM-like Higgs boson observed at the LHC. Here we find that the 
current limits are complementary, and give results that are very similar in reach, to those 
from direct detection experiments. 

The high mass regime, on the other hand, remains relatively unaffected by all the current 
bounds, and we find a viable region of parameter space for DM masses above 500 GeV. 
This provides an interesting scenario for direct detections experiments, since it was recently 
pointed out [H| that the future XENON- IT experiment should be able to probe essentially 
the entire viable IDM parameter space, including the high-mass regions that are likely to 
remain inaccessible to the LHC. 

Finally, we have studied the interplay between the dark matter phenomenology and the 
LHC Higgs phenomenology, in particular the possibility to have an enhanced h° 77 decay 
rate, as hinted in the present data. Here, the most important result of our analysis is that 
the latter can only be achieved for a dark matter relic density which is below the WMAP 
limit. Moreover, we find that the corresponding parameter points cannot be extrapolated to 
scales higher than 10 7 GeV with the requirement of perturbativity satisfied. Based on our 
findings, we have presented three benchmark scenarios that capture the essential features 
discussed throughout our analysis, i.e. dark matter relic density and detection, possibility 
to extrapolate the model to the GUT scale, and an enhanced h° — >• 77 decay rate inspired 
by the LHC discovery. 
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Appendix A: One-loop scalar masses 

In this Appendix, we give the expressions for the (one-loop) running scalar masses in the 
Inert Doublet Model, which have been obtained as discussed in Sec. [TTTl In the MS scheme, 
the masses for the SM-like Higgs boson (h°) and inert scalars (H° : A , are given by 



m 2 h o 



" W {* ( 6 + 3) MM2w) + ( 3 + 4) Ao(M|) 

+ c^M^ f 12 + B (M^, M 2 W ) + M*U + Bo (Ml M%) 



+ 



3m 



h° 



2M| 
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2M\) 



- 6Mf - 124^ 
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+ 2X 2 L v 2 B (m 2 H o,m 2 H o) + 2X 2 s v 2 B (m 2 A o,m 2 A o) + X 2 3 v 2 B (m 2 H ± , m 2 H ±) 



2„,2 D / 2 _2 



v2„,2; 



2 2 



a 



IQkc 2 w s 2 w 



A (M 2 ) + 2c 2 w A (M^) + 2c%7?% ± B (m H ± , Mj^) 



+ m AO B (m A o, M 2 ) - 2M 2 - 4c 2 w M^ 
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m 2 H ± = m 2 H ± 
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In these expressions is the Z boson mass, M w the W mass, and i/f (rrif) the Yukawa 
coupling (mass) of a generic fermion / that comes in number of colours (3 for quarks 
and 1 for leptons). We further denote the electromagnetic fine-structure constant as a ) 
and we use the shorthand notation = 1 — = M^/M z for the weak mixing angle. 
The functions A G and B are the regularized (finite) parts of the one- and two-point scalar 
integrals, respectively. The first argument p 2 = of the function B has been suppressed. 
For details on these integrals, as well as their expressions in terms of elementary functions, 
the reader is referred to Ref. 12811. 



Appendix B: One-loop beta functions 

The one-loop beta functions for the running of the parameters appearing in the general 
2HDM potential are known in the literature, see for example [311 ]. To collect everything 
using a consistent notation for the IDM, we give the results of our calculation here. The 
beta function for a running parameter \{(Q) is defined through 

167r2 ^g = ^ = ^ ) + ^ ) + ^ } ' (B1) 

where we choose to separate the contributions from scalars (s), gauge interactions (g) and 
Yukawa couplings (y). The contributions to the IDM beta functions from the scalar inter- 
actions are given by 

P^ = 24A 2 + 2A 2 + 2A3A4 + A 2 + A 2 , (B2) 
P% = 24A 2 + 2 A 2 + 2 A 3 A 4 + A 2 + A 2 , (B3) 
/3<f 3 } = 4 (Ax + A 2 ) (3A 3 + A 4 ) + 4A 2 + 2A 2 + 2A 2 , (B4) 
P { £ = 4A 4 (Ax + A 2 + 2A 3 + A 4 ) + 8A 2 , (B5) 
P® = 4A 5 (Ax + A 2 + 2A 3 + 3A 4 ) . (B6) 
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The contributions from gauge interactions read 



^ = I (3s 4 + g' 4 + 2<?V 2 ) - SXi (sg 2 + g' 2 ) , (B7) 

flU = I W + g' 4 + ig 2 g' 2 ) - ^ (s 9 2 + g' 2 ) , (bs) 

P$ = \ (3^ 4 + g' 4 - zg 2 g' 2 ) - sa 3 (s^ 2 + g' 2 ) , (bq) 

P$ = Zg 2 g' 2 - 3A 4 (3g 2 + g' 2 ) , (Bio) 

^ = -3A 5 (3/ + g' 2 ) . (Bll) 
Finally, including only the third generation fermions, we get from the Yukawa interactions 

< = 4Ax (y 2 + 3j£ + 3 % 2 ) - 2 (y< + 3 % 4 + 3yf) , (B12) 

^ = 0, (B13) 

f% = 2A 3 (y 2 T + 3y 2 + 3y 2 ) , (B14) 

^ = 2A 4 (y 2 + 3y 2 + 3y 2 b ) , (B15) 

= 2A 5 (y 2 + 3y 2 + 3y 2 ) . (B16) 

We also need the running of the gauge couplings, which is given by 

Pgt = b i9i (9i = 9, g', 9s), (B17) 

where, for the case of two Higgs doublets, b\ = —3, b 2 = 7, and 6 3 = —7. The corresponding 
expressions for the Yukawa couplings are 



= vt ( -y 2 g' 2 - \g 2 - 8g 2 + vl + \vl + \vl ) , (bis) 



f$ = y b (~g' 2 - \g 2 - + vl + \vl + \yl) , (bw) 

ffl = Vr ("^' 2 - -/ + \vl + 3% 2 + 3 % 2 ) . (B20) 
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